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ABSTRACT 



From the initial pure D2O phase, the SNO experiment has measured the flux of *B 
solar neutrinos via CC, NC, and ES interactions. The NC flux measurement, 0nc° ~ 
S.OQloiei X 10^cm~^s~^, represents the total flux of active neutrinos and is consistent 
with solar model expectations. When combined with the CC, 0cc*^ ~ l-'''6laio x 
lO'^cm-^s-i, and ES, = 2.39t[J:^^ x lO^cm-^s'^ results, SNO demonstrates the 
existence of neutrino flavour transformations at the 5.3cr level. 

The second salt phase of SNO operation is nearing completion. The de- 
tector calibration for this phase has been carried out and systematic uncertainties 
have been evaluated in preparation for the flrst salt publication. As demonstrated 
from MC simulation studies, 1 year of salt data livetime will produce CC and ES re- 
sults with similar precision as the pure D2O phase and a substantially improved NC 
measurement. In addition, precision flux measurements with no energy constraint 
will be possible. 
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1 Introduction 



The Sudbury Neutrino Observatory (SNO) [IJ, a 1000 T D2O Cerenkov detector lo- 
cated 2 km underground (6010 m water equivalent) in Sudbury, Canada, is sensitive 
to solar neutrinos via the interactions 

i^e + d ^p + p + e" (CC) 
i^x + d ^p + n + z/x (NC) 

i^ + e- ^ //x + e- (ES) (1) 

where i^x represents all active neutrino flavours. 

Whereas the charged current (CC) interaction is only sensitive to pq neutri- 
nos, the neutral current (NC) interaction is equally sensitive to all active neutrino 
flavours. This sensitivity to non-^e neutrinos, combined with the precise vq flux 
measurement from the CC channel, gives SNO the unique ability to test the hy- 
pothesis that a portion of the vqs produced in the sun undergo flavour change while 
in transit to earth. The elastic scattering (ES) interaction does have some sensitivity 
to non-i^e neutrinos but for SNO is a low statistics channel. 

Since the standard solar model predicts a negligible flux of non-z^e neutri- 
nos, a non-zero difference between the total flux of active neutrinos from the NC 
channel and the flux of electron neutrinos from the CC channel would provide a 
clear indication of neutrino flavour transformation. 

SNO has nearly completed the second of the flrst two phases of data taking. 
Results from the initial 'pure D2O' phase, in which the central target consisted solely 
of purifled heavy water, are described in the next section of this document. The 
flnal section provides an update of the status of the second 'salt phase' in which the 
characteristic response of NC events has been altered by the addition of NaCl to the 
D2O. 

2 Pure D2O Phase Results 

2.1 SNO Detector, Calibration, and Backgrounds 

The SNO Detector consists of 9456 photo-multiplier tube (PMT) and light concen- 
trator units arrayed on a geodesic support structure. This array surrounds a 1700 
T spherical shell of light water, which in turn encloses the spherical acrylic vessel 
containing the 1000 T of heavy water. 

Electrons produced in the D2O emit a cone of Cerenkov light which is 
detected by the PMT's. Neutrons produced in the detector can be captured by 
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a deuteron which subsequently emits a 6.25 MeV photon. Photons interact with 
electrons in the D2O via Compton scattering and these electrons are detected as 
described above. 

The event energy is approximately linearly related to the number of PMT's 
registering light. Vertex positions and event directions are reconstructed via fitting 
algorithms which use the relative PMT positions and firing times to reconstruct the 
Cerenkov cone. 

In order to understand the detector response and systematic uncertainties, 
a series of cahbration sources are placed at various positions within the detector vol- 
ume. Analyses of the calibration data and comparisons with Monte Carlo simulated 
events are then carried out. 

The optical properties of the detector are evaluated by transmitting laser 
light at wavelengths between 337-620 nm into the D2O volume via optical fibers 
affixed to a glass diffusing ball. The response of the PMT array, for various positions 
of the diffusing ball, is evaluated and the optical properties of the detector measured. 

In order to set the absolute energy scale and to evaluate energy scale and 
resolution uncertainties, an ^^N source|2], which emits 6.13 MeV photons, is utilized. 
These studies are supplemented by ®Li source [3] runs in which /9's with an end point 
near 14 MeV are produced and ^H(p,7)^He('pT') source |3] runs in which 19.8 MeV 
photons are emitted. In addition to the energy response, these sources are also used 
to determine the reconstructed event vertex and event direction detector response 
and uncertainties. 

The NC flux measurement requires an understanding of the detector neu- 
tron capture efficiency. A 252 source is introduced at various positions within the 
detector to produce a known quantity of neutrons and analysis of this data provides 
a measurement of the capture efficiency and uncertainties. 

Evaluations of backgrounds in the data set are carried out through a vari- 
ety of ex-situ and in-situ analyses. Instrumental backgrounds are evaluated through 
an examination of calibration data and a series of selection cuts adopted to sup- 
press them. Ex-situ measurements using ion exchange and membrane degassing 
techniques coupled with in-situ analyses of the SNO low-energy signal region are 
employed to estimate the level of radioactive thorium and uranium backgrounds. 
These backgrounds give rise to photons that contaminate the low-energy region 
of the SNO spectrum and, more crucially, can cause photo-disintegration of the 
deuterons inducing a neutron background signal indistinguishable from NC events. 

The combined systematic uncertainties from these studies are presented in 

Tabled 
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Source 


CC Uncert. % 


NC Uncert. 


% (t)^r Uncert. % 


Energy scale 


-4.2,+4.3 


-6.2,±6.1 


-10.4,±10.3 


Energy resolution 


-0.9,+0.0 


-0.0,±4.4 


-0.0,±6.8 


Vertex accuracy 


-2.8,+2.9 


±1.8 


±1.4 


D2O Cerenkov 


-0.1, +0.2 


-2.6,±1.2 


-3.7,±1.7 


PMT Cerenkov 


±0.1 


-2.1,±1.6 


-3.0,±2.2 


Neutron capture 


±0.0 


-4.0,±3.6 


-5.8,±5.2 


Total uncertainty 


-5.2,±5.2 


-8.5,±9.1 


-13.2,±14.1 



Table 1: The primary systematic uncertainties on the extracted charged current 
(CC) and neutral current (NC) neutrino fluxes are presented. In addition, the 
uncertainties on the measured non-electron neutrino flux {4>^t) are given. The total 
uncertainty for the elastic scattering flux (ES - not shown) is -4. 8, +5.0 percent. 

2.2 Pure D2O Phase Results 

After applying data cuts to reduce backgrounds, including a lower kinetic energy 
threshold cut of T> 5 MeV, and event vertex radius cut R< 550 cm, 2928 events 
are selected from the 306.4 days of live time collected between Nov. 2, 1999 and 
May 28, 2001. 

The reconstructed energy, angle between event direction and direction to 
the sun, and vertex radius spectra from this sample are simultaneously flt using an 
extended maximum likelihood technique to estimate the relative fractions of CC, 
NC, and ES events in the sample. 

The flt results, assuming a standard spectral shape jHj, are 1967.7lgJ;9 
CC events, 263.6l^|^ ES events, and 576.5l'^48;9 events; only statistical uncer- 
tainties are given. Figure Q shows the energy spectrum of the data overlayed by the 
flt distribution and the relative contributions from CC, ES, NC, and backgrounds. 

The measured fluxes of *B neutrinos from each of the three channels, in- 



cluding systematic uncertainties, are [ZI 

f^^^ = 1.76l°^(stat.)+S:S^ (syst.) X lOWV^ (2) 
<^|NO = 2.39+°|^(stat.)+S:;^ (syst.) x lO^cm-^g-i (3) 
^ 5.09^°;^^(stat.)+S:S (syst.) x lO^cm'^g-i^ 

Transforming these numbers into estimates of the flux of and the flux 
of u^T- (assuming only three active neutrino flavours) gives 

(Pe = 1.76l°:°^(stat.)+S:S^ (syst.) x lO^cm'^s'^ (5) 
(j)^r = 3.4ll°;^^(stat.)+!;:J (syst.) x lO^cm'^^-i^ 
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Figure 1: Shown are the pure D2O phase kinetic energy spectrum of neutrino events 
(points) and fit distribution (solid) (left) and confidence level regions in the 4>^T-'Pe 
plane (right). 

The measured non-zero flux of v^j^r neutrinos indicates, at a more than 5cr 
level, that vq produced in the sun undergo flavour transformation. Figure ^ indicates 
this result graphically, where it may be seen that the 0e measurement comes almost 
entirely from the CC result whereas the (p^r estimate is primarily given by the NC 
result. The ES result does provide additional information for both fluxes and is in 
good agreement with the CC and NC values. 

The NC measured total flux of active neutrinos given in Equation ^ 
is consistent with the standard solar model theoretical prediction jH] of (/'ssM ~ 
S.OSlgisi- However, this result requires constraining the shape of the energy 
spectrum. Relaxing this constraint produces the less precise result 

0SNO = QA2t\T,{stat.)t'o'i{sys.). (7) 

One of the leading theoretical candidates to explain neutrino flavour trans- 
formations is the MSW model jH]. SNO data has been evaluated in the context of 
this model producing the confidence level regions in the Am^-tan^ 6 plane given in 
Figure El P^- As can be seen, a large portion of MSW space is rejected at the 90% 
C.L. by SNO data. 

Combining the SNO measurements with results from CI [llj, Ga [T^ 1121 
IHj, and the Super-Kamiokande [12] experiments produces the confidence level con- 
tours given in Figure|21 As shown, current experimental evidence from solar neutrino 
experiments discounts all but LMA and LOW region solutions at better than a 99% 
C.L.. The relative minima in these two regions strongly favours the LMA solution. 
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= % ^min. =25.791 at (4.47e-01 , 4.73e-05) = 
10"' Vb8= 1.029 hep= 1.000 ^ 
= SNO pure D^O d/n spectra free = 

•^Q l ^ \ \ \ \ 

lo' 10^ 10^ 10"' 1 10 

tan^e 



l0""rX ^min.=69.917 at (4.07e-01 , 6.76e-05) ^ 
\ b8= 1 .035 hep= 1 .000 | 

10 'VSNO pure DjO d/n spectra -e 
- + SK-I zenith spectra + CI + Ga free : 

.|q"^^I \ \ \ \ \ 

lo" 10^ 10^ 10"' 1 10 

tan^e 



Figure 2: Shown are MSW C.L. region from SNO only pure D2O data (left) and for 
combined solar neutrino data (right). 



3 Salt Phase Update 

In May 2001, approximately 2 tonnes of NaCl was dissolved in the D2O volume of 
the SNO detector. The addition of the salt enhanced the neutron capture efficiency 
and the characteristic response to NC events. The first salt publication will focus 
on an initial data set corresponding to roughly 250 live days of running. 

3.1 Detector Calibration and Response 

As with the pure D2O phase, the detector response must be calibrated and sys- 
tematic uncertainties evaluated for the salt phase. Three examples highlighting the 
current status of calibrations in the salt phase are given here. 

The energy scale calibration is carried out via a combination of laser and 
^^N scans as outlined above. Full laser optical scans are taken at approximately 
3-4 month intervals. ^^N stability runs with the source positioned at the centre of 
the detector are taken every few weeks while full two-plane scans with the source 
located at many source positions throughout the D2O volume are taken every few 
months. 

Early in the salt phase it was discovered that the energy response of the 
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detector was changing with time. Figure El shows the mean effective number of 
PMT's (Neff) triggering per -"^^N event as a function of run date for -"^^N stabihty 
runs. As can be seen, the response changes by approximately ~ 2%/year. 

Subsequent evaluation of the laser data indicates that the measured D2O 
attenuations lengths were changing in a corresponding fashion. MC simulated 
^^N events have been generated assuming time varying attenuations and the mean 
effective number of PMT's for each MC run are shown in Figure 121 As can be seen, 
the change in response predicted by the MC for the measured changing attenuations 
agrees well with the measured change in the ^^N data. 

MC simulated electrons at a series energies are generated using the laser 
scan derived optical constants and a mapping between Neff and energy derived. 
After applying all corrections and this mapping function to the ^^N events, the 
mean energy response shown in Figure 01 is produced. Systematic uncertainties 
associated with the stability of detector response in time are evaluated by examining 
the variations in and differences between the data and MC ^^N energy distributions. 

Systematic uncertainties associated with spatial variations and asymme- 
tries in the detector are evaluated by examining ^^N data and MC runs taken 
at various positions throughout the D2O volume. Including additional uncertain- 
ties associated with source modeling, rate dependence, hardware stability, timing 
calibration, and cross-talk the total energy scale systematic uncertainty is roughly 
percent level. 

The primary method for distinguishing ES events from other signal and 
background events is to exploit the fact that the reconstructed event directions 
for such events will be correlated with the direction to the sun. Hence a good 
understanding of the detector angular response is required. 

Data and MC ^^N events are used to evaluate the systematic uncertainties 
on the angular response. What is important here is the agreement between data 
and MC spatially and temporally. 

For these events, the 'generated' event direction is taken to be the direction 
from the ^^N source position to the reconstructed event vertex. To ensure a reliable 
estimate for the 'generated' direction, a cut is placed requiring the reconstructed 
vertex of the electron to be at least 150 cm away from the ^^N source position. 

Figure |31 shows the distribution of the cosine of the angle between the fit 
and 'generated' direction for data and MC ^^N events for a run taken at the centre 
of the detector. In order to characterize the these distributions, the following simple 
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Corrected Nhit Mean vs SNO Julian Date 
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Figure 3: Mean effective number of PMT's triggering per ^^N event (left) for cen- 
trally positioned data (blue) and MC (red) stability runs and mean reconstructed 
energy (right) distributions as a function of run date are shown. 



three parameter function 

F(cos^) = Ar[e^2(i-cose) ^ p^gP3{i-cose)| 

is fit to the distributions. The overall normalization is arbitrary but left free in the 
fit. As can be discerned in Figure |2 there is good agreement between the data and 
MC fit parameters for this run. 

Figure |3] shows distributions of the three fit parameters for data and MC 
runs as a function of run date. Systematic uncertainties are derived from these distri- 
butions are dominated by the mean difference between the data and MC. Analogous 
distributions have been evaluated as a function of radial position within the D2O vol- 
ume with the total systematic uncertainty on these parameters assessed to be at the 
10-15% level. This translates into negligible uncertainties on the CC and NC mea- 
surements but does contribute substantially to the ES total systematic uncertainty. 
The ES measurement is, however, still dominated by statistical uncertainties. 

The most significant change with the inclusion of salt corresponds to the 
change in response to NC events (neutrons). In this configuration, neutrons are 
primarily captured on ^^Cl with a subsequent release of multiple photons 

n + =^^C7/ 36(7/ + S7 (Es-^ ~ 8.6Mey). (9) 

Thus not only does the higher capture cross-section on CI dramatically 



8 



Angle between Fit and 'Generated' Direction 



hadf 


hadf 


Entries 


11696 


Entries 


3086 


Mean 


0.7459 


l\1ean 


0.67 


RMS 


0.4586 


RMS 


0.5618 


/ ndf 


40.83 / 31 


x'indt 


29.59 / 31 


Prob 


0.1114 


Prob 


0.5387 


pO 


0.1829 ±0.005452 


pO 


0.1835 ±0.01 119 


pi 


0.2742 + 0.02628 


P1 


0.2796 ±0.04215 


P2 


22.65 ± 1.301 


p2 


25.8 ± 2.745 


p3 


4.795 ±0.2133 


p3 


4.922 ± 0.3493 



16 



Data - central 
MC - central ^^N 



N run 
run 



■ ni.i.i«k«,«,'i 



.1* 



1.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

cos(e) 





















, I \ ;1 




1 j 




















1 






; 4 ^ 















Figure 4: Shown are the distribution of the cosine of the angle between the fit and 
'generated' direction for data (blue) and MC (red) ^^N events for a run taken at the 
centre of the detector (left) and fit parameters for data (blue) and MC (red) vs run 
date. 



improve the neutron detection efficiency, but in addition the the multiple photon 
nature of these events causes them to exhibit a greater degree of isotropy than, for 
example, CC events with a single Cerenkov cone. 

Figure El compares CC and NC MC simulated events for isotropy parameter 
9ij. This is the mean pair angle between PMT's that have triggered in an event. As 
can be seen, such information is clearly useful in separating CC and NC events. In 
fact, utilizing this information will allow a precise determination of the fluxes to be 
made without the use of the energy spectrum. 

The neutron efficiency has been measured via analysis of neutrons produced 
from a 252Qf source placed at various positions within the D2O volume. Figure 
shows the comparison of neutron detection efficiencies between the pure D2O and 
salt phases as a function of radial position within the D2O. The response at the 
centre of the detector is almost a factor of 2 higher in salt but, owing to the greater 
difference at high radii, the detection efficiency inside R=550 cm is improved by 
almost a factor of 3. 



3.2 Signal Extraction 

Signal extraction for the salt phase is carried out in the same fashion as pure D2O 
phase. The neutrino candidate event distributions are fit, using an unbinned ex- 
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Figure 5: Shown are comparisons between CC (blue) and NC (red) MC 9ij isotropy 
distributions (left) and between pure D2O (red) and salt (blue) phase neutron detec- 
tion efficiency distributions vs radial position (right). 



tended maximum likelihood technique, with NC, CC, and ES signal pdf 's plus with 
fixed background pdf's included. The primary difference for the salt phase is the 
altered shape of the NC pdf's and the inclusion of an additional isotropy variable. 

In order to evaluate the expected statistical sensitivities of the salt phase 
dataset, signal MC simulated samples with rates corresponding to the measured 
pure D2O fluxes were generated flt with various combinations of salt pdf's. One live 
year of data was simulated and cuts corresponding to the pure D2O selection were 
applied. 

Table El presents the percent statistical error from the published pure D2O 
results and for the salt MC study. When all four variables (energy, radius, angle 
between flt direction and direction to sun, and isotropy) are included in the flt, the 
CC and ES uncertainties are projected to be similar to the pure D2O but the NC 
measurement improves by almost a factor of 2. 

Interestingly, when the energy pdf's are not included in the flt the CC 
uncertainty does not increase appreciably and the NC uncertainty remains substan- 
tially better than in the pure D2O phase. Hence with the salt phase data a precision 
measurement will be made with no assumptions on the shapes of the signal energy 
distributions. 
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Table 2: Statistical uncertainty MC projections for salt phase. Variables are 
E— energy, R— event radius, Ogun—o-ngle between fit direction and direction to sun, 
and (3u—isotropy parameter. 



Phase 


Variables in Fit 


CC Error 


NC Error 


ES Error 


Pure D2O 




3.4% 


8.6% 


10% 


Salt MC 




4.2% 


6.3% 


10% 


Salt MC 


E,R,^sun,/3l4 


3.3% 


4.6% 


10% 


Salt MC 




3.8% 


5.3% 


10% 



4 Summary 

Precise measurements of the CC and NC fluxes from the pure D2O phase of SNO 
running have been presented. These results can be translated into values for electron 
neutrino and non-electron neutrino fluxes. The measurement of a statistically signif- 
icant non-electron neutrino flux provides strong evidence for flavour transformation 
of electron neutrinos produced in the sun and can be consistent with an MSW oscil- 
lation model. The shape constrained measurement of the total *B neutrino flux is 
consistent with standard solar model predictions. The unconstrained measurement, 
while still consistent, is dominated by statistical uncertainties. 

Various aspects of detector calibration and response for the salt phase have 
been presented including energy scale, angular resolution, and neutron response. As 
indicated, the salt phase detector response is well understood and these analyses are 
in an advanced stage in preparation for the forthcoming initial salt phase publication. 

MC studies indicate that 1 year of salt livetime will produce a CC result 
with similar statistical precision as the pure D2O phase and an NC result with 
substantially improved statistical precision. Note that the currently pubhshed NC 
measurement uncertainty is dominated by statistics. In addition, it has been demon- 
strated that a precision energy-unconstrained flt can be made with the salt phase 
data. 
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